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ABSTRACT 

We present an analysis of the redshift-space power spectrum, P{k), of rich clusters 
of galaxies based on an automated cluster catalogue selected from the APM Galaxy 
Survey. We find that P{k) can be approximated by a power law, P{k) oc fc", with 
n K, —1.6 over the wavenumber range 0.04 h Mpc"^ < fc < 0.1 h Mpc"^ Over this 
range of wavenumbers, the APM cluster power spectrum has the same shape as the 
power spectra measured for optical and IRAS galaxies. This is consistent with a simple 
linear bias model in which different tracers have the same power spectrum as that of 
the mass distribution but shifted in amplitude by a constant biasing factor. On larger 
scales, the power spectrum of APM clusters flattens and appears to turn over on a scale 
k ~ 0.03 h Mpc^^. We compare the power spectra estimated from simulated APM 
cluster catalogues to those estimated directly from cubical N-body simulation volumes 
and find that the APM cluster survey should give reliable estimates of the true power 
spectrum at wavenumbers k > 0.02 h Mpc~^. These results suggest that the observed 
turn-over in the power spectrum may be a real feature of the cluster distribution and 
that we have detected the transition to a near scale-invariant power spectrum implied 
by observations of anisotropics in the cosmic microwave background radiation. The scale 
of the turn-over in the cluster power spectrum is in good agreement with the scale of 
the turn-over observed in the power spectrum of APM galaxies. 
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1 INTRODUCTION 

Most modern theories of structure formation, including 
those based on inflation or topological defects in the early 
Universe, predict a near scale-invariant spectrum of den- 
sity perturbations on large spatial scales. On small spatial 
scales, A < 10 h'^Mpc, the observed fluctuations in the 
galaxy distribution difi^er very substantially from a scale- 
invariant form. For example, the power spectra estimated 
from infra-red and optically selected redshift surveys are rea- 
sonably well approximated by P{k) oc k~^'^ at wavenumbers 
k > 0.1 h Mpc~^ and there is no convincing evidence for a 
turn-over to a scale-invariant form, P{k) oc fc, at smaller 
wavenumbers. However, such a turn-over must exist since 
measurements of the microwave background fluctuations on 
angular scales > 1° suggest a near scale-invariant spectrum 
of perturbations on spatial scales > 100 h~^Mpc (see the 
reviews by Scott, Silk & White 1995 and Bond 1996). 

A convincing detection of a turn-over in the power spec- 
trum of density irregularities is therefore of fundamental sig- 
nificance and would establish a direct link between fluctua- 
tions observed in the microwave background radiation and 



those observed in the density distribution. Furthermore, in 
theories such as the Cold Dark Matter (CDM) model, a peak 
in the power spectrum is predicted on the scale of the Hub- 
ble radius at the time that matter and radiation have equal 
density, X^qu ~ 10(n/i^)"^Mpc (see e.g. Efstathiou 1990). 
An observation of a peak in the power spectrum of galaxy 
clustering can therefore be used to constrain the parameter 
r — Qh that is fundamental to CDM models. 

There is some tentative evidence for a turn-over at 
a wavenumber k ~ 0.025 in the three-dimensional power 
spectrum inferred from the two-dimensional clustering of 
galaxies measured from the APM galaxy survey (Baugh & 
Efstathiou 1993, Baugh & Efstathiou 1994, see for exam- 
ple Figure 11 of Baugh & Efstathiou 1994). Gaztanaga & 
Baugh (1997) have carefully investigated the signiflcance 
of this observed turn-over by repeating, using simulations 
of the APM galaxy survey, the procedure of inverting the 
two-dimensional clustering statistics to obtain the three- 
dimensional power spectrum . They conclude that the turn- 
over, which occurs between 0.02 < fc < 0.06 h Mpc~^, can 
be reproduced by the inversion process. However, at these 
small wavenumbers, the systematic errors in the construe- 
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tion of the APM Galaxy Survey are difficult to quantify and 
could perhaps be greater than the random errors (see Mad- 
dox, Efstathiou & Sutherland f 996 for a detailed discussion 
of systematic errors in the APM survey), ft is therefore dif- 
ficult to assign a meaningful statistical significance to this 
result. 

The detection of a peak in the power spectrum of the 
galaxy distribution is one of the main scientific goals of the 
Anglo- Australian 2dF redshift survey (see e.g. Efstathiou 
1996) and the Sloan Digital Sky survey (Gunn & Weinberg 
1995) which aim to measure redshifts of ~ 10® galaxies. 
However, in this paper we show that redshift surveys of 
much smaller numbers (~ 10^) of rich clusters of galaxies 
can provide accurate measurements of the power spectrum 
on large spatial scales. Relatively little work has been done 
on the three-dimensional power spectrum of rich clusters of 
galaxies (see Gramman & Einasto 1992, Peacock & West 
1992, Einasto et al. 1993). Most recent work has focused on 
measurements of the two-point spatial correlation function 
of rich clusters of galaxies, ^(r), and whether systematic er- 
rors in the cluster catalogues affect the amplitude of ^(r) 
on scales < 20 /i~^Mpc. The closest investigation to the 
one presented here is that of Peacock and Nicholson (1991), 
who describe a power spectrum analysis of a redshift survey 
of 310 radio galaxies. Their results are described in further 
detail in Section 4. 

In this paper we analyse the redshift survey of 364 APM 
clusters described by Dalton et al. (1994). In previous papers 
(Dalton et al. 1992, Dalton et al. 1994), we have applied a 
number of tests to show that the APM cluster catalogue is 
free of the projection and selection biases known to affect 
clustering in the Abell cluster catalogues (Sutherland 1988, 
Efstathiou et al. 1992). As we show here, the large volume 
surveyed by the APM cluster survey (V > 3 X lO'^/i-^Mpc^) 
renders it suitable for an investigation of clustering on very 
large scales. 

This paper is laid out as follows. In Section 2 we summa- 
rize the techniques used in the power spectrum analysis. We 
apply these techniques to the APM cluster survey and inves- 
tigate the sensitivity of our results to the selection function, 
weighting function and background cosmological model. In 
Section 3 we construct simulated APM cluster catalogues 
from large N-body simulations to quantify any systematic 
errors and biases in our power spectrum estimator. From 
this analysis we can assess the significance of the observed 
peak in power spectrum of APM clusters. Our conclusions 
are presented in Section 4. 



2 ESTIMATION OF THE POWER SPECTRUM 

The analysis described here is similar to the power spec- 
trum analysis of the Stromlo-APM galaxy redshift survey 
described by Tadros & Efstathiou (1996) (hereafter TE96). 
We follow their notation unless otherwise stated. We ap- 
ply the power spectrum analysis to the redshift survey of 
APM clusters described by Dalton et al. (1994). We use the 
largest statistically uniform sample of APM clusters. Sam- 
ple B of Dalton et al. (1994) containing 364 clusters over 
the southern APM area (21'' < RA < S**, -72.5° < dec < 
— 17.5°). The number density of clusters in this sample is 
~ 3.4 X 10-^( h-^Mpc)-^. The redshift distribution of the 
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Figure 1. Histogram of the redshift distribution of the 364 clus- 
ters in Sample B of the APM cluster redshift survey. The lines 
shows the result of smoothing the histogram with a Gaussian of 
width u for four values in the range 3000 - 10000 kms~^ as 
indicated in the figure. 



sample is shown in Figure ^. The median redshift of the 
sample is Zmed ~ 0.09 and the cluster distribution extends 
to a redshift z ~ 0.2. 

2.1 Method 

We treat the APM cluster redshift survey in the same way 
as a flux limited galaxy survey, and use the methods of Feld- 
man. Kaiser & Peacock (1994) (hereafter FKP), as imple- 
mented by TE96, to account for the radially varying selec- 
tion function. To define the selection function of the survey 
we smooth the velocity distribution shown in Figure ^ with 
a Gaussian of width v kms~^. The fiducial value of v is 
4000 kms"\ 

The estimation of the power spectrum may be summa- 
rized as follows (for further details the reader is referred to 
TE96). We compute a weighted density field 



F(r) 



w (r) [uc (r) — aus (r)] 



) (r) d'^rl 



(1) 



where the subscript c denotes the cluster density in the real 
catalogue and s denotes the density field for a random cata- 
logue with same angular and radial selection functions as the 
cluster survey. In equation (1), n(r) is the expected mean 
density of clusters in a catalogue with the same angular and 
radial selection functions as the data. The radial selection 
function is derived from the fits to the redshift distribution 
plotted in Figure ^. The function n (r) can be separated into 
the mean galaxy density n (r) as a function of radial distance 
r, multiplied by the angular mask of the catalogue {i.e. the 
distribution of the APM galaxies on the sky defined by the 
survey boundary and excluding regions around bright stars, 
globular clusters etc.). The factor a is the ratio of the space 
densities in the real catalogue to that in the random cata- 
logue. In the analysis presented here we use several thousand 
times as many points in the random catalogue as there are 
clusters in the real catalogue, and compute a from the ratio 
of the sums 



E 



(1 -f 47rn(n)J3) 



(2) 



where we have set 47rJ3 = 40000(/i ^Mpc)'^ and the sums 
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Figure 2. Three-dimensional power spectra of the APM clusters redshift survey. Panel (a) shows the effect of using different values of 
Pw{^) in the weighting function of equation (3). Three values are used: P„(fc) = 50000, 92000 and 130000 ( /i-^Mpc) as shown by the 
crosses, circles, and squares respectively. The estimates in Panel (a) assume a spatially flat geometry with = 0. Panel (b) shows the 
effect of varying the assumed cosmology. The power spectrum is shown for two spatially flat models: Qrn = 1-0, Qa = 0.0 (circles), and 

f2m = 0.2, Ca = 0.8 (triangles). Both spectra in the right hand panel have Pw{k) = 92000 ( h~^yipcj in the weighting function. For 
clarity, sets of points have been offset from each other by one seventh of the bin spacing in both panels. 



run over all clusters and random points (see Efstathiou 1996, 
Section 5.3 for details). The specific choice for J3 is based 
on the power spectra of CDM models multiplied by a factor 
~ 4 to account for the bias of clusters with respect to the 
underlying mass distribution. However, the results presented 
below are insensitive to the exact value of AnJa used in the 
analysis. 

Each cluster in the analysis is weighted according to the 
minimum variance weight function (FKP) 



which depends on the value of Pw(k) at each wavenumber 
k. We have set (fc) in equation (3) equal to a constant 
value and investigated how the estimates of the power spec- 
trum change for three values of Pw (k): 50000, 92000 and 
130,000 ( ft-^Mpc) , which span the range of the observed 
power spectrum over a wide range of wavenumbers. The 
power spectrum is derived by Fourier transforming equation 
(1), removing shot noise (equation 2.1.9 of FKP) and aver- 
aging the power-spectrum estimates over shells in fc-space of 
volume Vk ■ The power spectrum estimates will be correlated 
over a range in fc-space given approximately by 5k ~ 
where D is the characteristic depth of the survey. To obtain 
roughly independent points, all power spectra presented in 
this paper (except those calculated from the full periodic 
box of the N-body simulations described in Section 3) have 
been binned over a range in wavenumber ~ . 



To perform the power spectrum analysis one must de- 
fine a background cosmological model. Firstly, a cosmology 
must be assumed in the conversion of cluster redshifts to 
proper distances, and secondly a cosmological volume ele- 
ment must be used to obtain the radial number density from 
the solid line in Figure |^. Since the median redshift of clus- 
ters in the APM survey Zmed = 0.09, the assumed cosmology 
may affect the results. Consequently we have carried out the 
analysis for two assumed cosmologies. Firstly a spatially flat 
universe with a matter density parameter fi^ = 1-0, and 
secondly a spatially flat universe with fim = 0.2, and a cos- 
mological constant contributing most of the closure density, 
JIa = 0.8. The results for an open universe with S7m = 0.2 
and JIa = are intermediate between those for the two spa- 
tially flat models. 



2.2 Computation of Errors 

Error bars on P(/c) are computed using equation (2.4.6) of 
Feldman, Kaiser & Peacock (1994). Their expression was 
derived assuming that the density field of the tracer objects 
followed a Gaussian distribution. This will not be true on 
small scales where the density fluctuations in the cluster 
distribution are greater than unity, even if the primordial 
fluctuations were strictly Gaussian. Thus, in Section 3.2 we 
test the accuracy of the error estimates using mock cluster 
catalogues generated from large A^-body simulations. 
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2.3 Results for P{k) 

In Figure ^ we show the power spectrum for the 364 clus- 
ters in the APM cluster redshift survey. The cluster sample 
was truncated at a proper distance of rmax ~ 600 /i~^Mpc 
and the selection function for the sample was derived by 
smoothing the cluster redshift distribution with a Gaussian 
of width 4000km s~^ as illustrated in Figure 1. 

Figure 2(a) shows the cluster power spectra cal- 
culated assuming values of Pw{k) = 50000,92000 and 
130, 000 ( /i^^Mpc) in the weighting function (equation 
3) and distances computed assuming a spatially flat uni- 
verse with zero cosmological constant. Figure 2(b) shows 
the effect of changing the underlying cosmological model 
whilst keeping the weighting parameter fixed at Pwik) = 
92, 000 ( /i~^Mpc) . The circles show the power spectrum 
derived assuming a spatially flat cosmology with fim = 1.0, 
Q.A = and the triangles show the power spectrum calcu- 
lated assuming fim = 0.2, SI a = 0.8. This comparison shows 
that uncertainties in the background cosmological model can 
affect the amplitude of the cluster power spectrum by a fac- 
tor of up to ~ 1.3, but that the shape of the power spectrum 
is not significantly affected. Figure 2 shows that the power 
spectrum of APM clusters is well approximated by a power 
law P(k) (X fc" with index n = —1.6 ±0.3 over the wavenum- 
ber range 0.04 < fc < 0.1 h Mpc~^. The power spectrum 
flattens at smaller wavenumbers and appears to turn-over a 
wavenumber k ~ 0.03 h Mpc~^. 



2.4 Sensitivity of P(fc) to weighting and 
selection function 

Figure ^ shows that the power spectrum of APM clusters 
is not very sensitive to the value of Pw{k) assumed in the 
weighting function of equation (3). Changing the weighting 
factor Pu, (fc) by relatively large amounts leads to changes in 
the power spectrum points that are very much smaller than 
the error bars. 

The power spectrum estimates are more sensitive to the 
to the value of u, the velocity smoothing parameter that is 
used to generate a smooth selection function from the cluster 
redshift distribution (see Figure 1). Figure 3(a) shows power 
spectra calculated using three values for this smoothing, 
3000km s~^ and 5000km s^^ (spanning the fiducial value of 
4000km s~^ used in the estimates of Figure 2) and a much 
larger value, v = 10000km s~^. For these estimates we use 
a constant weighting factor of Pn>{k) = 92000 ( /i~^Mpc) . 
Although the power spectra are insensitive to the smooth- 
ing velocity for small changes around the fiducial value of 
4000km s"\ if we smooth with v = 10, 000km s ^ the am- 
plitude of the power spectrum at wavenumbers k IS 0.035 
rises sufficiently to eliminate the turnover in P{k). For the 
largest smoothing, the selection function is a poor fit to ob- 
served redshift distribution, and particularly to the small 
number of high redshift clusters {z > 0.13, see Figure |l|). 
We find that the sensitivity of the power spectrum esti- 
mates to the selection function smoothing can be reduced 
substantially by truncating the sample at a proper distance 
of Vmax = 400 /i~^Mpc, thus eliminating the high redshift 
tail of the cluster n{z) distribution. This is illustrated in 
Figure 3(b) which shows P{k) estimates for the same three 
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Figure 3. Figure 3(a) shows the effect on the power spec- 
trum estimates of changing the value of u used to smooth the 
cluster redshift histogram (see Figurel). Estimates of P{k) are 
shown for three values of u: u = 3000 kms"'^ (crosses), u = 
5000 km s~^ (circles) and v = 10000 kms~^ (squares). The 
power spectrum becomes sensitive to the value of u for wavenum- 
bers k < 0.035 h Mpc-i. Pi gure 3(b) is identical to Figure 
3(a) except that we have limited the cluster sample to clus- 
ters with proper distances of less than 400 h~^Mpc. For all of 
these estimates, we have assumed a constant weighting factor of 
Pw{k) = 92000 ^ h^^Mpcj and a spatially flat background cos- 
mology with A = 0. For clarity, sets of points have been offset 
from each other by one seventh of the bin spacing in both panels. 



velocity smoothings as in Figure 3(a) but where we have 
truncated the cluster sample at rmax = 400 /i~^Mpc. In 
this figure the three estimates of P{k) are compatible to well 
within the la error bars. The efi'ect of the velocity smoothing 
is to produce a slight change in the amplitude of P{k) and 
all three estimates show a turnover at fc ~ 0.03 h Mpc~^. 

Unless otherwise stated, in the rest of this paper we 
will use the APM cluster power spectrum computed under 
following assumptions: 

(i) P{k) = 92000 ( h'^MpcY in the weighting function of 
equation (3). 

(ii) A smoothing velocity u = 4000 km s~^ to define the 
selection function. 

(in) r^ax = 400 h'^Mpc. 

(iv) A spatially flat background cosmology with zero cosmo- 
logical constant. 

To summarize, this section shows that the shape of the es- 
timated power spectrum is insensitive to these assumptions 
to within the la error bars. However, plausible variations 
of these assumptions can lead to a change in the overall 
amplitude of P{k) by up to a factor of 1.5. 
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3 COMPARISON WITH SIMULATED APM 
CLUSTER SURVEYS 

In this section we investigate the reliability of the power 
spectrum estimator described in Section 2.1. We measure 
P{k) from mock APM cluster surveys extracted from large 
N-body simulations to determine the minimum wavenumber 
at which the true power spectrum can be recovered accu- 
rately. From this investigation, we can assess the significance 
of the turn-over seen in the power spectrum of APM clusters 
at a wavenumber of ~ 0.03 h Mpc~^ (Figures 2 and 3). 



3.1 Construction of simulated surveys 

We have constructed mock APM cluster redshift surveys 
from A'^-body simulations of two cold dark matter (CDM) 
models. The methods for creating mock catalogue from 
the numerical simulations are described in more detail by 
Croft and Efstathiou (1994). The simulations consist of 
two ensembles of 9 simulations each containing 160^ par- 
ticles within a periodic computational box of length lb = 
600 Ti-^Mpc. They are similar to the simulations described 
in detail by Croft and Efstathiou (1994) but employ more 
particles within a larger computation box. (The simulations 
of Croft and Efstathiou use 100'^ particles within a box of 
size lb = 300 /i~^Mpc). The simulations were run with the 
particle-particle-particle-mesh (P"^M) code described by Ef- 
stathiou et al. (1985) and model gravitational clustering 
in a CDM dominated universe with scale invariant initial 
density fluctuations. The two ensembles are as follows: the 
standard CDM model (Davis et al. 1985), i.e. a spatially 
flat universe with fio = 1 and h = 0.5 (the SCDM en- 
semble); a spatially flat low density CDM universe with 
Jlo = 0.2, h — 1.0, and a cosmological constant contribution 
= A/(3Ho) = (1 - ^o) = 0.8 (the LCDM ensemble). 
The initial power spectra of the models are generated from 
the fitting function 



P(k) oc 



1 -I- (^ak + (bk) '^ + [cky 



(4) 



where u = 1.13, a = (6.4/r) h'^Mpc, b = (3.0/r) h'^Mpc 
and c = (1.7/r) /i~^Mpc. Equation (4) is a good approxi- 
mation to the linear power spectrum of scale-invariant CDM 
models with low baryon density, fit <^ (Bond & Efs- 
tathiou 1984). The parameter F in equation (4) is equal to 
Qoh. Thus F = 0.5 for the SCDM ensemble and F = 0.2 for 
the LCDM ensemble. 

The final output times of the models are chosen to ap- 
proximately match the microwave background anisotropics 
measured in the first year COBE maps (Smoot et al. 1992) 
ignoring any contribution from gravitational waves. Thus 
the rms mass fluctuations in spheres of radius 8 h~^Mpc 
are as = 1. The clustering statistics of clusters of galaxies 
are a very weak function of time, so the simulation output 
time used is not a critical consideration (see e.g. Croft & 
Efstathiou 1994). 

Candidate cluster centres are selected by using a 
friends-of-friends algorithm to link particles with separation 
less than 0.1 times the mean inter-particle separation. These 
group centres are then used as the centres of spheres of size 
Vc — 0.5 /i~^Mpc within which the mass of the cluster is 



calculated. The centre of mass of the cluster was found and 
any cluster which was less than rc from the centre of an- 
other, richer cluster was deleted. The procedure was then 
repeated, and a final list of clusters was generated and or- 
dered by mass. A lower mass cut off was applied to generate 
a cluster sample of a chosen mean number density. 

The mock surveys extracted from the cluster distribu- 
tion are designed to sample the same volume as the APM 
cluster survey using the same sky mask as the real data {i.e. 
the identical area of the sky and excluded areas such as holes 
around bright stars and globular clusters). The lower mass 
cut-off of clusters identified in the N-body simulations was 
adjusted to match the space-density of clusters in the real 
APM cluster redshift survey. The mock catalogues were lim- 
ited at a proper distance of 600 /i~^Mpc; however we have 
checked that the results and conclusions of this section are 
unchanged if the mock catalogues are limited at the smaller 
distance of 400 /i^^Mpc. The smoothed APM cluster selec- 
tion function {ly — 4000km s~^) was imposed on the cluster 
catalogues and (since the peculiar velocities of the clusters 
in the N-body simulations are known) we shifted the clus- 
ters to their redshift space positions to create mock APM 
catalogues in redshift space. 



3.2 P{k) from simulated surveys 

Figure ^ shows the power spectra measured from the mock 
APM cluster catalogues, using the estimator described in 
Section 2. The crosses show the average value of P{k) mea- 
sured from 27 mock surveys for each ensemble using three 
randomly chosen centres for each simulation. The error bars 
show the standard deviation of the mean P{k), computed 
from the scatter over the 27 mock surveys. The circles 
in Figure ^ show the mean redshift-space power spectrum 
for each ensemble computed from the full A'^-body simula- 
tion volume {i.e. with no application of a selection func- 
tion or angular sky mask). The values of Pw{k) used in 
the weighting function were Pw{k) = 50000 ( /i~^Mpc) and 

P„(fc) = 16000 ( h-^MpcY for the LCDM and SCDM sim- 
ulations respectively. These values of Pw{k) were chosen to 
provide a match to the mean amplitude of P{k) from these 
two cosmological models over the range of wavenumbers 
plotted in the figure. The figure shows that our estimator 
of P{k) based on equation (1) provides a good measure of 
the power spectrum of APM-like cluster catalogues in the 
LCDM model for wavenumbers k > 0.02 h Mpc"^ For the 
SCDM model, we see some differences in the power spectra 
estimates at wavenumbers 0.02 - 0.05 h Mpc~^, but these 
are not large enough to eliminate the break in the power 
spectrum estimates for the mock APM catalogues. 

TE96 give a detailed analysis of biases in estimating 
the power spectrum from redshift surveys. There are two 
effects that can lead to systematic differences between the 
estimated and true power spectra on large scales: (a) the es- 
timated power spectrum is a convolution of the true power 
spectrum and the power spectrum of the window function 
of the survey; (b) the estimated power spectrum is biased 
low at small wavenumbers if the mean space density of the 
sample is estimated from the survey itself. Both effects de- 
pend on the geometry of the survey and on the form of the 
true power spectrum and will affect the estimates of P{k) 
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Figure 4. The circles plotted in each of the two panels show the redshift-space power spectrum of clusters in the LCDM model (left hand 
panel) and the SCDM model (right hand panel) computed from the full periodic volume of the N-body simulations. The crosses show the 
power spectra estimated from the average over 27 realizations of the APM cluster redshift survey constructed from the simulations. The 
error bars in all cases show the error on the mean value of P{k). The power spectrum estimated from mock APM cluster catalogues is in 
good agreement with the power spectrum computed from the full simulation volumes on scales k > 0.02 h Mpc~^ for the LCDM model, 
and k > 0.03 h Mpc~^ for the SCDM model. The solid lines show the linear theory spectra for the two models normalized arbitrarily to 
match the data points at wavenumbers k ~ 0.05 h Mpc~^. The inset in the left hand panel shows the performance of the error estimator 
(equation 2.4.6 of FKP). The triangles in the inset plot show the error bar on the mean P{k) obtained from the 27 realizations of the 
mock APM catalogues. The crosses show the error, o"(A:), derived by using equation (2.4.6) of FKP. The error bar plotted on each of the 
crosses shows the size of the error on the mean error, obtained from the spread in the measured cr{k) over the 27 realizations of the APM 
cluster catalogue. 



from the mock APM catalogues plotted in Figure y. How- 
ever, these biases are not present in the power spectrum 
estimates from the full periodic volumes of the N-body sim- 
ulations (plotted as the circles in Figure ^). Thus, the close 
agreement between the power spectra in each of Figures ^(a) 
and ^ (b) shows that neither bias is significant for wavenum- 
bers k > 0.02 h Mpc"^ 

The inset in the left hand panel of Figure ^ shows a 
test of the error estimator used in this analysis. The crosses 
show the error bars calculated from equation (2.4.6) of FKP, 
averaged over 27 realizations. The errors on each of these 
points show the error on the mean. The triangles show the 
error on P{k) calculated by evaluating the spread in the 
measured value of P{k) over the 27 realizations of the APM 
cluster survey. Thus equation (2.4.6) of FKP provides an 
excellent estimate of the true error on P{k) for all but the 
largest wavenumbers plotted in the figure. 

In Figure ^ we compare the power spectrum of the APM 
clusters with the power spectra measured from the mock 
catalogues described above. The points for the real survey 
differ slightly in the two panels because in each case we as- 
sume background cosmologies consistent with the LCDM 
and SCDM theoretical models (c/ Figure 2b). The APM 
cluster power spectrum cannot be matched by the SCDM 
model, in agreement with our previous analyses of the spa- 
tial two-point correlation function of the APM cluster sur- 



vey (Dalton et al. 1992, Dalton et al. 1994). The LCDM 
model is much closer to the observations in shape and am- 
plitude, but under-predicts the clustering strength and does 
not reproduce the turnover in the observed power spectrum 
at fe < 0.03 h Mpc~^. A scale-invariant CDM model with 
a lower value of F (F ~ 0.15) would likely provide a bet- 
ter match to the clustering amplitude of rich clusters at 
wavenumbers k ~ 0.04 h Mpc~^ (e.g. see Figures 8 and 
9 of Gaztanaga, Croft & Dalton 1995) but would fail even 
more strongly to match the downturn in the observed power 
spectrum at small wavenumbers. 

These discrepancies between theory and observation do 
not seem particularly problematic, given the uncertainties 
in both the observational estimates and theoretical mod- 
els. Neither the SCDM or LCDM models considered here 
provide an acceptable match to the observations, but these 
belong to a highly restricted class of CDM model. It may be 
possible to achieve a better match with a CDM-type model 
by allowing other parameters to vary, e.g. by incorporating 
a tilt in the primordial spectral index (White et al. 1995), 
or by invoking an admixture of hot and cold dark matter 
(Klypin et al. 1993). 
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Figure 5. A compaj-ison of the power spectrum of APM clusters with the predictions of the LCDM and SCDM models. The crosses 
in each case show the average power spectra of 27 mock APM cluster surveys extracted from the N-body simulations, together with Icr 
errors appropriate to a single mock survey. The circles show the power spectra computed for the real APM cluster survey assuming a 
background cosmology that is consistent with the theoretical models. 
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Figure 6. The power spectrum of APM clusters together with that of IRAS galaxies, optical galaxies and radio galaxies. The galaxy 
power spectra have been scaled upward by a scale-independent factor (see Table 1) to match the amplitude of the APM cluster power 
spectrum. 
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Table I 





k h Mpc"^ 




be 




P 


IRAS 


0.05 - 


0.1 


6.43 


± 0.73 


1.65 


0.18 


Stromlo-APM 


0.05 - 


0.1 


3.13 


± 0.32 


0.19 


0.90 


Radio gals 


0.04 - 


0.1 


0.99 


± 0.15 


1.34 


0.25 



Table 1. Values for the relative bias (squared) between IRAS, optical and radio galaxies compared to APM clusters. Column 2 gives 
tile range in wavenumber over which the fit was performed. Column 3 gives the factor by which the galaxy power spectrum must be 
scaled to match the clusters power spectrum i.e. the relative bias squared. The last two columns give the P^i' degree of freedom and 
the probability for the fit to a scale-independent relative bias. 



4 DISCUSSION 

We first compare the power spectrum of rich clusters with 
those for other types of object in the Universe. Figure ^ 
shows the APM cluster power spectrum together with power 
spectra for optical Stromlo-APM galaxies (Loveday et al. 
1992a, TE96), IRAS galaxies (from the combined QDOT 
and 1.2Jy surveys, as calculated in Tadros & Efstathiou 
1995) and a sample of 310 radio galaxies (as calculated by 
Peacock & Nicholson 1991). The spectra for optical, radio 
and IRAS galaxies have been scaled upwards by constant 
factors bl to match the amplitude of the APM clusters power 
spectrum. The factors be are thus the relative bias parame- 
ters for the three samples and are computed by performing 
a weighted fit to the ratio of the observed power spectra over 
a specified range in wavenumber. Table pj lists the range in 
wavenumber over which the data were fitted, the values of 
the relative bias factors be, as well as the probability for 
the fit. A reasonable value for the fit probability shows that 
these power spectra are consistent with the hypothesis of a 
linear relative bias between the various tracer objects over 
the scales indicated in Table 1. 

Peacock & Dodds (1994) found a relative bias param- 
eter of 2.37 between radio galaxies and Abell clusters. This 
differs from the lower relative bias between radio galaxies 
and APM clusters listed in Table 1. This is a result of the 
higher amplitude of the power spectrum of Abell compared 
to APM clusters. However, there is now a large body of 
evidence to show that the clustering of Abell clusters is 
enhanced by non-uniformities in the Abell catalogue and 
does not reflect the true clustering properties of rich clus- 
ters (Sutherland 1988, Efstathiou et al. 1992, Croft et al. 
1997). The relative bias of be — 0.99 between radio galaxies 
and APM clusters from Table |l| is thus consistent with the 
idea that luminous radio galaxies preferentially inhabit the 
cores of rich clusters of galaxies. 

The agreement between the shapes of the power spec- 
tra for rich clusters of galaxies, IRAS and optically selected 
galaxies is consistent with a simple linear bias model. This 
provides an argument against models in which the galaxy 
formation process introduces spatial correlations so that 
galaxies are non-linearly biased with respect to the mass 
distribution (see e.g. Babul & White 1991). Figure ^ sug- 
gests a simpler interpretation, in which different types of 
galaxies and rich clusters are linearly biased on large-scales 
and each traces the shape of the underlying matter power 
spectrum. 

We have found some evidence from the power spectrum 
of the APM cluster sample for a peak at fc ~ 0.03 h Mpc~^, 




0.01 0.1 



k (h Mpc 

Figure 7. The filled circles show the power spectrum of APM 
clusters as plotted in Figure 6. The open squares show the 
power spectrum of APM galaxies estimated by inverting the two- 
dimension angular correlation function as described by Maddox 
et al. 1996 (see their Figures 32 — 35). The APM galaxy power 
spectrum has been multiplied by a factor of 5 to match the am- 
plitude of the cluster power spectrum. The error bars were deter- 
mined from the scatter in the inversions from four nearly equal 
area zones of the APM survey. The two lines show linear theory 
power spectra for a scale-invariant LCDM model with Ha = 0.8, 
h = 1 (solid line) and a mixed dark matter model with fli, = 0.3, 
h = 0.5. The linear spectra have been normalised arbitrarily to 
match the cluster power spectrum approximately in the region 
k ~ 0.05 h Mpc-i. 



followed by a downturn at smaller wavenumbers. The tests 
described in Section 3 suggest that the biases in the esti- 
mator of P{k) are too small to cause such a downturn. Fig- 
ure 7 shows the APM cluster power spectrum (filled circles) 
together with the three-dimensional power spectrum deter- 
mined by inverting the angular two-point correlation func- 
tion measured for the APM Galaxy Survey (open squares). 
The latter points are from the analysis of Maddox et al. 
(1996) who applied the inversion technique developed by 
Baugh and Efstathiou (1993). The amplitude of the APM 
galaxy power spectrum has been multiplied by a factor of 
5 to match the amplitude of the cluster power spectrum. 
Both the cluster and galaxy power spectra show a peak at a 
wavenumber ~ 0.03-0.04 h Mpc~^, which suggests that the 
peaks reflect a real feature of the underlying mass distribu- 
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tion. As mentioned in the introduction, such a peak must ex- 
ist at a wavenumber > 0.01 h Mpc~^ since the anisotropies 
of the cosmic microwave background radiation on angular 
scales > 1° are well approximated by a scale-invariant ini- 
tial spectrum (Bond 1996). 

The lines in Figure 7 show linear theory power spectra 
for two models, normalized to match the observations at a 
wavenumber k ~ 0.05 h Mpc"\ The V = 0.2 LCDM model 
has too much power on large scales to provide a good match 
to eithe r the cluster or galaxy power spectra. The dashed line 



shows . . scale-invariant mixed dark matter (MDM) power 



spectrum, with parameters = 0.3, Q.h = 0.05, — 0.65, 
and h = 0.5 (where il,^, Qi, and Qrn are the respective cos- 
mological densities in massive neutrinos, baryons and cold 
dark matter). The power spectrum of the MDM model is 
computed from the fits given in Ma (1996), but over the 
range of wavenumbers plotted in Figure 7 there is little dif- 
ference between the MDM spectrum and a CDM spectrum 
of the form (4) with F = Qmh ~ 0.3. As Figure 7 shows, 
a CDM-like model with F ~ 0.3 has a peak in the power 
spectrum that approximately matches the observations. 

Because of the low space density of rich clusters of 
galaxies and their enhanced clustering strength compared 
to normal galaxies, rich clusters provide a powerful probe 
of large-scale structure in the Universe. At wavenumbers 
k > 0.01 h Mpc"\ redshift surveys of only a few thousand 
rich clusters are capable of providing comparable results on 
the power spectrum, to those from redshift surveys of order 
10^ galaxies. Figure 7 thus suggests that a modest increase 
in the size of the cluster redshift sample could confirm the 
reality of a peak in the power spectrum at a high significance 
level. 

Acknowledgments: We thank John Peacock for com- 
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